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We calculate the energy spectrum of an electron moving in a two-dimensional lattice which is 
defined by an electric potential and an applied perpendicular magnetic field modulated by a periodic 
surface magnetization. The spatial direction of this magnetization introduces complex phases into 
the Fourier coefficients of the magnetic field. We investigate the effect of the relative phases between 
electric and magnetic modulation on band width and internal structure of the Landau levels. 

Energy spectra of electrons moving in the presence of a perpendicular homogeneous magnetic field B through a 
two-dimensional (2D) lattice are determined by fascinating commensurability effects. For a 2D lattice described by 
\Q ■ the electrostatic potential 

V(x, y) = V x cos(Kx) + V y cos{Ky) (1) 
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with lattice constant a — 2n/K, the single-particle energy spectrum can be calculated if the magnetic flux <E> = Ba 
through a unit cell of the lattice is a rational multiple of the flux quantum <j> = ch/e, 

$/$o ee 2ir/K 2 l 2 =p/q. (2) 



Here I = (ch/eB) 1 / 2 is the magnetic length, and q and p are relative prime integers. If the amplitude of V(x,y) is 
' much smaller than the distance between Landau levels, i.e. the cyclotron energy ?wj c ee heB/mc, mixing of Landau 
levels can be neglected. Then the resulting energy spectrum consists of broadened Landau bands, each with the same 
internal structure which, when plotted versus the inverse flux ratio &o/Q, is known as Hofstadter's butterfly Q], see 
Fig. 1(a). The width of these Landau bands oscillates as a function of the Landau quantum number n and the flux 
ratio $/$o- 

So far the gap structure predicted by Hofstadter's butterfly has not been observed directly in experiments, but 
there is strong evidence that it plays an important role for the magnetoresistance of the laterally modulated 2D 
q | electron systems in GaAs-AlGaAs heterostructures (for a review see Ref. ||). Magnetoresistance measurements have 
been performed on ID and 2D lateral superlattices, with characteristically different results. In ID superlattice so 
called "Weiss oscillations" are observed, which can be understood from the oscillating width of the Landau bands. 
The dispersion of these bands leads to a group velocity, the quantum analog of the classical guiding center drift of 
cyclotron orbits in the oscillating electric field. As a consequence, an oscillatory contribution to the conductivity 
results, which is called "band conductivity" and which is absent in homogeneous 2D electron systems. In a ID 
superlatttice the Weiss oscillations of the band conductivity can also be calculated quasi-classically from the 
guiding center drift. In weak 2D superlattices this band conductivity is drastically suppressed, and this is attributed 
to the splitting of the Landau bands into narrow subbands separated by large gaps, as predicted by the Hofstadter 
spectrum. In such a spectrum the group velocities and, as a consequence, the band conductivity are much smaller 
than in the corresponding ID case. 

The analog of the Weiss oscillations in 2D electron systems with a periodic magnetic field, which had been predicted 
H and investigated in some detail several years ago, has been experimentally observed [p|—[To|| only recently. The 
magnetic modulation was achieved by metallic strips on the surface, consisting of superconducting || or ferromagnetic 
P,|To[ material and exerting at low temperatures stress on the underlying semiconductor, so that the magnetic mod- 
ulation was found to be accompanied by a distinct electric modulation. Systematic investigations [ pl]Jl~2| have shown 
that this stress-induced electrostatic modulation always has a considerable content of higher harmonics. Nevertheless, 
the experimental results for unidirectional magneto-electric modulation are in reasonable agreement with the theo- 
retical calculations j(j,|7j which, in the regime of the Weiss oscillations may even be replaced by their quasi-classical 
limit 

Quite recently, magnetoresistance in 2D magneto-electric superlattices has been investigated Q,^fj . The lattice was 
defined by a quadratic array of ferromagnetic posts deposited on the surface of a GaAs-AlGaAs heterostructure, so that 
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a combined magnetic and electric modulation results at low temperatures. The observed oscillations show distinct 
features which can not be explained by the quasi-classical theory, which is sufficient to explain the corresponding 
results for ID modulations. Therefore the question arises if, similar to the case with purely electric modulation, a 
Hofstadter-type energy spectrum suppresses the modulation-induced band conductivity and how Hofstadter's butterfly 
is modified. We will investigate this question for the most simple modulation model which includes the experimentally 
relevant symmetries. 

In the experiment, the ferromagnetic posts can be magnetized parallel or perpendicular to the plane of the 2D 
electron gas. As a simple model for the magnetic field resulting for an arbitrary direction of the periodic magnetization, 
we consider the field of identical magnetic dipoles m located at the lattice sites (n x a, n y a, 0) in the plane z = 0. It is 
a straightforward exercise to calculate the Fourier expansion of the resulting periodic magnetic field in a plane z < 0, 
where the 2D electron gas is assumed to be located. For the z component one obtains 

B z {r;z) = Y,Bt{z)e^\ (3) 

9 



where 



B^z) = — (im-q + m z q) e~ q ^ , (4) 



with q = \q\, q = (n x K,n y K,0) and integers n x and n y . Thus the complex phase of the Fourier coefficients B^(z) 
is determined by the orientation of the magnetic moments on the surface, a result which remains valid beyond this 
simple dipole approximation. Also the exponential damping of higher Fourier components is typical for magneto- and 
electrostatic fields satisfying Laplace's equation. 

We now assume a strictly 2D model for the electron gas and neglect the in-plane components of the periodic 
magnetic field together with the z dependence of B z (r). To describe B z (r) in the plane of the electrons, we take 
Eq. (0) for the fixed z value at this plane and use in the 2D Hamiltonian the vector potential A(r) with Fourier 
coefficients Aq — i{q y , — q Xl 0)B^/q 2 , which commutes with the momentum operator. For the homogeneous external 
magnetic field B we use the Landau gauge A^ (r) = (0, xBq, 0), so that, in the presence of an electrostatic potential, 
the Hamiltonian can be written as 

H = H {0) + H m + H cm , (5) 

where 



ff(o) = J_ 
2m 



Pi + (Py + -xBq] 



(6) 



is the Hamiltonian of the homogeneous system with the usual Landau eigenenergies En = huj c (n+l/2) and eigenstates 
|n, k y ) centered around xq = ~l 2 k y . The other two terms are Fourier sums as in Eq.(||). For the term linear in the 
magnetic modulation we have 

= -g [q y p x ~ q x {p y + muj c x)] ioq, (7) 

where we have introduced u)q = eB^/mc in analogy to the cyclotron frequency uj c = eBo /mc of the homogeneous 
system. The quadratic term in the magnetic modulation can be rewritten as a single Fourier expansion, and then has 
the same form as the elctrostatic modulation, to which it adds as 



H% m = V s + M s , (8) 



with 



m (q — k) ■ k ,„. 

M ff = > — — ^ t u t , (9) 

q 2 ^ ( g --fc) 2 fc2 «- & k ' y > 

where the sum is over all k — (m x K,m y K,0) with integers m x and m y . Note that u>^—o = whereas M^ =0 yields a 
constant energy shift which we neglect. The matrix elements of H ora in the Landau basis are immediately given by 
the known ^] matrix elements (n', k' y \ exp(i<f • r)\n, k y ). Those of H m can also be reduced to these matrix elements 
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after the action of the momentum and position operators in Eq. (|7|) on the Landau functions has been evaluated using 
the recurrance relations of the Hermite polynomials. It is easily seen that the Hamiltonian of Eq.(|j) has the same 
symmetry under lattice translations and gauge transformations as the Hamiltonian in the case with purely electric 
modulation. As a consequence, in order to calculate its eigenstates and -energies, we may take the same steps as 
in the treatment of the purely electric Hofstadter problem 0]. First, we assume both the electric and the magnetic 
modulation to be weak enough, so that the coupling of different Landau levels may be neglected. Then, within one 
Landau level, we make for the eigenstates the ansatz 

oo 

||ra,a}) = 22 c x( n ,Oi)\n,ky + XK) (10) 

A— — oo 

with the restriction \k y \ < K/2. The coefficients c\(n,a) must be determined from an eigenvalue equation 

oo 

{H x ,,x-E n!a 5 x >,\}c x {n,a) = 0. (11) 

A— — oo 

The matrix elements of H are of the form H$\ — 3\',\En \ and, for the modulation-induced terms, 

K,x = E *A',a+£ S n(Q) e- il2 ^ k v +XK+S ^ (12) 
<? 

with 

S™(q) = G n (Q)hwj, (13) 



Sr(q)=£n(Q)(V ? +Mf). (14) 

Here we use the notation [§ Q = l 2 q 2 /2, C n (Q) = exp(-Q/2) L n (Q), and Q n {Q) = exp(-Q/2) [L { £_ x {Q) + L n {Q)/2], 

where Ln\Q) is an associated Laguerre polynomial [L n = Ln \ = 0]. Equation jll| ) is a generalization of 

Harper's equation, which contains arbitrary Fourier components of the periodic magnetic and electrostatic fields. 

We now retain only the simplest Fourier contributions necessary to describe typical experimental situations. First, 
we assume that the (stress induced) electrostatic potential has inversion symmetry with respect to the center of a 
magnetized post, i.e. has a pure cosine expansion with real coefficients Vx. Second, we consider only a simply harmonic 
modulation of the magnetic field in both directions, u)^ x = lo x for g x — (K, 0) and uj = tu y for g y = (0, K), and, of 
course = ut, and lo^ = for all other values of q. This implies finite contributions for with q = ±2g x and 
±2<7 a , and Mq — for all other values of q. Finally, we retain in the electrostatic potential energy only contributions 
with the same wavevectors, i.e. the ground and the second harmonics V±g lt and V±2$ for fi — x or y, and we assume 
that these coefficients are real. Then, with E 7ha — hu> c (n + 1/2) + e„ iQ , the generalization ( |Tl] ) of Harper's equation 
reads 



2 Re (S n (g x ) e^ K l k y+ XK i + S« m (2g x ) e -^K[k y +XK^ 
+S n (g y )c x+1 + S n (g y )*c x -i + S^(2g y )c x+2 + $T (2&)* CA _ 2 - , (15) 

where 

S n (g») = V^C n {Q) + hcjf, g n (Q) , (16) 

for fi = x or y, and Q = l 2 K 2 /2. 

Let us first ignore the second-harmonic terms S^ l m (2g tl ). Then Eq. ( |j~5"| ) takes the usual form of Harper's equation 
however with complex coefficients iSn(<?/i) = l^l exp(z<y9 M ). If one introduces c\ = cxp(—i\ip y )c\, the c\ satisfy 
Harper's equation with the real coefficients \S y \ in the off-diagonal and 2\S X \ cos(l 2 K\k y + XK] — ip x ) in the diagonal 
terms. 

Thus, for a modulation in x direction only, i.e. for Sy = 0, one obtains simple cosine bands of width 4 1 tS^ | . This 
generalizes results obtained by Peeters and Vasilopoulo s M for the special case of (i) real uj x , i.e. a cosine expansion of 
the magnetic modulation field ( "in-phase modulation" M) as can be realized according to Eq. (^) by a magnetization 
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in z direction, and (ii) imaginary u> x , i.e. a pure sinus expansion of B z {x) ( "out-of-phase modulation" Q) as can be 
realized by a magnetization parallel to the plane of the 2D electron system. 

For a square-symmetric modulation, S n {g y ) = S n (g x ), the same bandwidth factor iS^I as for the ID modulation 
can be factorized out of Harper's equation. If the commensurability condition (j^) holds, a Bloch ansatz reduces the 
infinite set ( |l5| ) to a closed set of p equations. The energy eigenvalues can be written in the form [|| 

en,a = \S X \ e(k x , k y ;j) , (17) 

where e(k x , k y ;j) consists of p subbands (J=l,...,p) defined in the first magnetic Brillouin zone \k x \ < ir/qa, \k y \ < 7r/a 
and is the same for all Landau bands n, except a possible shift in the repeated zone scheme due to the n-dependent 
phase ip x . The spectrum of the internal energy structure e(k x ,k y ;j) plotted versus the inverse flux is thus just 
Hofstadter's butterfly, Fig. 1(a). For a purely magnetic modulation this result was already known ||l6). We see that, 
in the purely harmonic approximation, the relative phase between magnetic and electric modulation determines only 
the width of the Landau bands, which is proportional to the n-dependent modulus of S x . The internal structure is 
not affected by this phase. 

We now discusss the effect of the second- harmonic terms in Eq. ([IB]) . According to Eq. (j|) the magnetic modulation 
always introduces such terms, which may, however, be small for a weak, purely harmonic magnetic modulation. But 
a second-harmonic electric modulation, which is known to be important for stress-induced modulations, will have the 
same effect. A nontrivial phase difference between first and second order terms may result, for instance, if the second 
order is dominated by an electric modulation, which we assumed to have real Fourier coefficients, while the first order 
is dominated by a magnetic modulation due to an in-plane surface magnetization. Therefore, it is of principle interest 
to investigate how these terms may change the butterfly. If we write <S° m (2<7 M ) = {S^l exp(2ix /i ) and express the 
c\ in Eq. ( |i~5| ) in terms of the c A , the coefficient of c A+2 becomes \S y 2\ exp(2zAx y ) and the corresponding diagonal 
term becomes 2 |«Sa;2 1 cos 2(l 2 K[k y + \K] — <p x ~ A\ x ), with Ax p = x M — cp^. In Fig. 1 we have plotted some typical 
spectra for the case of square symmetry, S y 2 = S X 2, where only a single relevant phase difference Ax remains. The 
normalization of the spectra is the same as in Fig. 1(a), i.e. the amplitude of the second-harmonic terms is chosen as 
\Sxi/S x \ = Cexp(— 37r$o/2 ( I > ), and C is taken constant, appropriate for the lowest Landau level, n — 0. Figure 1(b) is 
for C — 0.5 and A^ = 0. The exponential factor suppresses the second-harmonic term as the inverse flux ratio $o/$ 
approaches unity, the spectrum reduces to that of the unperturbed butterfly in Fig. 1(a). This factor also destroys 
the reflection symmetry about the value &o/& = 0.5. If we replace A\ = by A\ = 7r/2, i.e. change the sign in 
all second-harmonic terms, the spectrum is reflected about the line of zero energy (measured with respect to hto c /2). 
The spectrum in Fig. 1(c) is for A\ = tt/2, but for a large magnitude of the second-harmonic term, C = 2, and the 
asymmetry about the energy zero is enhanced. Note the strong crossing of subbands for c E>o/ < I > < 0.5, which destroys 
the large gap in the lower left part of Fig. 1(a). Figure 1(d) is for C — 2 and A\ = 7r/4. It is nearly, but not exactly, 
reflection symmetric about the zero energy line, as we found at this phase difference also for other values of C. Thus 
we see that the relative phase between the Fourier coefficients of the fundamental modulation and its second harmonic 
does not only affect the (n-dependent) widths of the Landau bands, but also manifestly their internal structures. 

Although the details of the energy spectrum, and as a consequence of the density of states, are changed by the 
second-harmonic terms, and in general their symmetry is lowered, the overall gap structure remains similar. Thus, 
qualitative consequences of the gap structure, for example the suppression of the band conductivity in 2D superlattices 
as compared with ID modulations, should remain valid in superlattices defined by mixed electrostatic and magnetic 
modulations even in the presence of second harmonics. Other harmonics and the coupling of different Landau levels, 
which have been investigated for purely electrostatic modulations [^],[l7| and have the tendency to lower the symmetry 
of the Hofstadter-type spectra, are not expected to change the picture qualitatively. 
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acknowledges financial support by the Deutsche Forschungsgemeinschaft. 
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FIG. 1. Scaled energy spectra, e n ,a/S x , versus inverse flux ratio $o/$ for a mixed magnetic and electric modulation 
with square symmetry, in the n — Landau level. The relative strength of the second-harmonic terms is given by 
\S X 2/S X \ = Cexp(— 37r$o/2<I>), the relative phase by Ax = \ axg(S X 2/S x ). The parameter values are (a) C = 0, Ax = 0, 
(b) C = 0.5, A X = 0, (c) C = 2.0, A X = tt/2, (d) C = 2.0, A X = tt/4. 
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